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Compressively strained SrTa0 2 N thin films were epitaxially grown on SrTi0 3 substrates using nitrogen 
plasma-assisted pulsed laser deposition. Piezoresponse force microscopy measurements revealed small 
domains (lO^lO 2 nm) that exhibited classical ferroelectricity, a behaviour not previously observed in 
perovskite oxynitrides. The surrounding matrix region exhibited relaxor ferroelectric-like behaviour, with 
remanent polarisation invoked by domain poling. First-principles calculations suggested that the small 
domains and the surrounding matrix had tons- type and a cis- type anion arrangements, respectively. These 
experiments demonstrate the promise of tailoring the functionality of perovskite oxynitrides by modifying 
the anion arrangements by using epitaxial strain. 



Perovskite oxynitrides have attracted much attention because of their novel electronic functionalities, such as 
visible light absorption 1 , photocatalytic activity 2 , colossal magneto-resistance 3 , and high dielectric constant 4 . 
These features presumably originate from the concomitant interaction between oxygen and nitrogen ions, 
which have different ionicities/covalencies 5 . The electronic functionalities of perovskite oxynitrides may also be 
influenced by the geometrical configuration of the O and N atoms around their metal cations. For A J B0 2 N-type 
oxynitrides, in which each B cation is surrounded by four O and two N ions, there are two possible anion 
configurations: the two nitrogen ions can occupy either adjacent (ds-type) or opposite (trans-type) sites in a 
£0 4 N 2 octahedron. Researchers have argued that the dielectric properties of A50 2 N are related to anion 
arrangement 6 " 8 : for example, Page et at. suggested that ferroelectricity in trans-type anion-ordered ATa0 2 N 
(A = Sr and Ba) phases may be caused by the off-centre displacement of Ta ions. They investigated this con- 
cept by theoretically studying the stability of phases with different nitrogen arrangements and space groups 6 . 
However, trans-type phases in this system are less energetically stable than ds-type phases 6 , and bulk SrTa0 2 N 
specimens have been confirmed to exhibit ds-type configurations (Fig. la) 7,8 . As suggested by Page et at. 6 , 
the ferroelectric trans-type phase would be stabilised by compressive epitaxial stress because of induced tet- 
ragonal distortion (Fig. lb). Such metastable ferroelectricity was not observed in an epitaxial thin film of 
compressively strained BaTa0 2 N deposited on a SrRu0 3 -buffered SrTi0 3 (STO) substrate 9 . This lack of ferroe- 
lectricity was likely caused by the considerably relaxed BaTa0 2 N lattice, indicated by the small lattice distortion 
(c/a = 1.0099) generated by the large lattice mismatch (—4.7%) between the film (a = 4.1 125 A 4 ) and the SrRu0 3 
buffer layer (a = 3.92 A; pseudo-cubic approximation 10 ). To apply enough tetragonal distortion to produce 
ferroelectric behaviour, SrTa0 2 N is a more suitable film than BaTa0 2 N, because the lattice constant of SrTa0 2 N 
(a = 4.03 A; pseudo-cubic approximation 11 ) produces a better lattice match ( — 3.1%) to the STO (a = 3.905 A) 
substrate. 

In the present paper, we report room-temperature ferroelectricity in compressively strained SrTa0 2 N thin 
films epitaxially grown on STO single- crystal substrates. Local piezoresponse measurements revealed the incorp- 
oration of small classical ferroelectric domains with dimensions of lO^lO 2 nm in a relaxor ferroelectric-like 
matrix. The ferroelectric domains and the relaxor-like matrix were identified as trans- and ds-type phases, 
respectively, based on results of first-principles density functional theory (DFT) calculations that mapped the 
total energy of SrTa0 2 N for various lattice parameters. 
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Figure 1 | Schematic illustration of anion arrangement in SrTa0 2 N. 

(a) Unstrained bulk SrTa0 2 N with a ds-type anion arrangement 7 ' 8 . 

(b) Tetragonally strained SrTa0 2 N thin film with a trans-type anion 
arrangement epitaxially stabilised on the SrTi0 3 substrate. The structures 
were drawn by VESTA 27 . 

Results 

Epitaxial growth of SrTa0 2 N thin films. SrTa0 3 _ x N x thin films 
were grown on the (100) planes of STO single-crystal substrates 
using nitrogen plasma- assisted pulsed laser deposition (NPA- 
PLD). The nitrogen content, x, in these films, was evaluated by 
nuclear reaction analysis (NRA); x had a range from 0 to —1.2, 
depending on the partial pressure of the supplied nitrogen gas and 
the deposition rate of the film (Supplementary Fig. SI). X-ray 
diffraction (XRD) confirmed that pure perovskite SrTa0 3 _ x N x 
films grew epitaxially on the STO substrates despite the wide 
variation in nitrogen content (Supplementary Fig. S2). The 
tetragonal lattice parameters, a and c, and the unit cell volume, V, 
of the SrTa0 3 _ x N x thin films monotonically increased with 
increasing x (Fig. 2a and 2b). Near the stoichiometric composition 
[x = 1.02(6)], we observed a V of 65.385 A 3 , agreeing well with the 
reported lattice volume of bulk SrTa0 2 N (65.478 A 3 ) 4 . The cation 
stoichiometry of this film was confirmed by Rutherford backscat- 
tering spectrometry (RBS) (Supplementary Fig. S3). In addition, we 
confirmed the band gap energy of 2.27 eV (Fig. 2c) determined from 
a plot of (ochv) 2 versus hv, which agreed well with a previously 
reported value (2.3 eV) 12 . Later, we shall discuss the crystal 
structure and physical properties of stoichiometric SrTa0 3 _ JC N JC 
thin films with nitrogen compositions of x = 1.00 ± 0.06, which 
we refer to hereafter as SrTa0 2 N. 

Investigating the microstructure of the SrTa0 2 N films using trans- 
mission electron microscopy (TEM; Fig. 2d) and high- angle annular 
dark field scanning TEM (HAADF-STEM; Fig. 2e), we found a sharp 
film/substrate interface and a uniform perovskite lattice structure with 
no segregation. Using XRD, we investigated the tetragonal lattice dis- 
tortions of the SrTa0 2 N thin films. Despite the partially relaxed lattice 
of the thin films (the representative reciprocal space map shown in the 
inset of Fig. 2b), the tetragonal distortion ratio of cla reached 1.026 {a 
= 3.98 A; c = 4.08 A) at a film thickness of 20 nm. The cla ratio 
decreased with increasing film thickness, becoming nearly constant 
(—1.013) in the films thicker than 100 nm. These results indicate that 
the SrTa0 2 N epitaxial thin films were highly strained and distorted, 
especially near the film/substrate interface, compared with bulk 
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Figure 2 | Structure of the SrTa0 3 _ x N x epitaxial thin films, (a) In-plane 
(circles) and out-of-plane (triangles) lattice constants of the SrTa0 3 _ x N x 
thin films as a function of nitrogen content, x. The solid lines are guides to 
the eye, representing a linear relationship, (b) Lattice volume calculated 
from (a). The dashed line shows the lattice volume of bulk SrTa0 2 N 4 . The 
inset shows an XRD reciprocal space map of the 103 diffraction of the 
SrTa0 2 N thin film and STO substrate. The dotted lines represent the lattice 
constant of bulk SrTa0 2 N, derived from Ref. 4. (c) Absorption coefficients 
calculated from the extinction coefficients. The left inset is a plot of (ochv) 2 
versus hv (solid line) and linear fitting results (dashed line). The right inset 
is a photograph of a patterned SrTa0 2 N thin film grown on a SrTi0 3 
substrate demonstrating visible light absorption, (d) Low-magnification 
cross-sectional TEM image of the SrTa0 2 N film, (e) High-magnification 
HAADF-STEM image of the SrTa0 2 N film. 

SrTa0 2 N {cla - 1.002) and BaTa0 2 N epitaxial thin films {cla - 
1.009) 9 . We also confirmed that SrTa0 2 N epitaxial thin films were 
perfectly relaxed {cla — 1.002) on (LaAlO 3 ) 0 .3(Sr 2 AlTaO 6 ) 0 .7 substrate 
{a = 3.868 A) with larger lattice mismatch (-4.0%) than STO. 

Investigation of ferroelectricity in SrTa0 2 N epitaxial thin films. 

The ferroelectricity of the SrTa0 2 N thin films was investigated using 
contact-resonance mode vertical piezoresponse force microscopy 
(VPFM) 13 , which probes the out-of-plane piezoresponse of the 
film. VPFM measurements were conducted under vacuum 
(<10 Pa) at room temperature. In the virgin state, small domains 
(size range from tens to hundreds of nanometers) in the thin film 
produced much larger piezoresponse than did the surrounding 
regions (Fig. 3a and 3b). To investigate these locally varying 
ferroelectric properties, we performed switching spectroscopic 
measurements, in which we recorded the piezoresponse of the 
remanent polarisation as a function of the amplitude of the preced- 
ing switching DC voltage pulse (Fig. 3c) 14 . The curves of local 
piezoresponse versus DC bias exhibited clear hysteresis (Fig. 3d) in 
the domains that exhibited large piezoresponse in the virgin state 
(point A in Fig. 3b). The amplitude components of the curve were 
butterfly shaped (Fig. 3e) and the phase components exhibited sharp 
180° inversions at the coercive fields (Fig. 3f). These features denote 
the presence of classical ferroelectricity, which has never before been 
observed in perovskite oxynitrides to our knowledge. The surround- 
ing regions with negligible piezoresponse in the virgin state (point B 
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Figure 3 | Local ferroelectric behaviour of the SrTa0 2 N epitaxial thin 
film, (a) Topographic and (b) piezoresponse (A cos0) images of the 
SrTa0 2 N thin film (thickness = 288 nm). (c) Schematic illustration of the 
switching DC pulses and probing AC waves in switching spectroscopy 
measurements. After each DC pulse, the piezoresponse originating from 
the remanent polarisation is probed with a small AC bias, (d)-(i) Local 
piezoresponse examined by switching spectroscopy with an AC bias of 
±1 V at positions A [(d)-(f)] andB [(g)-(i)] indicated in (b). Open circles 
represent the initial point of the measurements, (e), (f) and (h), (i) present 
the amplitude and phase components of (d), (g), respectively. 

in Fig. 3b) exhibited an apparent hysteretic piezoresponse induced by 
the external DC bias (Fig. 3g-3i). The piezoresponse of the SrTa0 2 N 
films was much weaker than those of conventional ferroelectric 



materials. The weaker response may have been caused by the small 
piezoelectric constant of SrTa0 2 N, but it may also be possible that 
the region near the film/substrate interface was only ferroelectric 
because of lattice relaxation, as discussed in the previous section. 

To investigate the unusual hysteretic behaviour observed in the 
matrix region, we took a series of VPFM images from the same region 
over time after performing domain poling 1516 . In the virgin state 
before poling, the piezoresponse signals were hardly detectable 
(Fig. 4b) because the measured area was much larger than the size 
of the ferroelectric domains (tens to hundreds of nanometers). 
Immediately after poling, the induced piezoresponse was clearly dif- 
ferent between the positively and negatively poled regions (Fig. 4e), 
and there was a 180° phase difference between the two distinct 
domains, indicating perfect polarisation inversion (Fig. 4c and 4d). 
The induced polar domains began to relax within 1 h, and the con- 
trast between the domains weakened considerably after 8 h (Fig. 4f- 
4i); in contrast, the surface morphology did not change over time 
(Fig. 4a). Such relaxation behaviour is typical in relaxor materials and 
has been observed in lead lanthanum zirconate titanate ceramics 17 . 

Interpreting these results must be done with caution, however: 
recent PFM studies on oxides reported that even paraelectric materi- 
als, such as Ti0 2 thin films and heterostructures such as LaAl0 3 / 
SrTi0 3 , can exhibit spurious relaxor-like behaviour, caused by reman- 
ent piezoresponses generated from external- electric field- induced 
motion of ionic species 18,19 . To check for this possibility, we examined 
how the dielectric constants of our thin films depended on tempe- 
rature (Fig. 5). We found no significant signs of ferroelectric- 
paraelectric phase transitions up to 600 K, at which point the sample 
began to be oxidised, as confirmed by XRD (data not shown). 
However, the real part of the dielectric constant, e', from 20 Hz to 
10 kHz, exhibited positive de'/dT, as seen in conventional "intrinsic" 
relaxors 20 . In contrast, in "extrinsic" relaxors, de'/dT tends to decrease 
with increasing temperature 21 . Thus, we tentatively conclude that the 
majority of the present SrTa0 2 N film surrounding the small domains 
(i.e., the relaxor-like matrix region) is an intrinsic relaxor. 

To investigate the possibility of ferroelectric-paraelectric phase 
transitions in the matrix region, we also conducted high-temperature 
XRD. To avoid oxidisation, the SrTa0 2 N film was capped with 
—20 nm of amorphous LaAl0 3 . Figure 6 shows the measured lattice 
constants as a function of temperature, clearly indicating an inflec- 
tion at —600 K in both the a- and c-axis lengths. Choi et al. reported 
similar behaviour in BaTi0 3 ferroelectric thin films at their Curie 
temperature 22 , leading us to speculate that the inflection point in the 
present SrTa0 2 N epitaxial thin film corresponds to its Curie 
temperature. 

Based on these findings, we believe that our SrTa0 2 N film con- 
tained classical ferroelectric domains (lO^lO 2 nm) and a surround- 
ing matrix with a relaxor-like nature, although we observed 
hysteretic behaviour only in the microscopic PFM measurements 
and could not obtain more direct evidence for ferroelectricity from 
macroscopic P-E measurements because of the large leakage current 
under high DC electric field (data not shown). 

Discussion 

First, we discuss the origin of the classical ferroelectricity in the 
present films. It is plausible that the domains exhibiting classical 
ferroelectricity can be classified as the polar trans-type phase, con- 
sidering that room-temperature ferroelectricity in the ds-type struc- 
tures has never been observed, experimentally or theoretically. 
However, first-principles calculations performed by Hinuma et al. 
suggested that SrTa0 2 N with specific ds-type structures exhibits 
reversible polarisation, although the polarisation flips within picose- 
conds at room temperature 23 . Such a fluctuating polarisation might 
be frozen by the epitaxial strain from the STO substrate, but our ab 
initio molecular dynamics simulations excluded the possibility of 
polarisation freezing (Supplementary Fig. S4). Thus, we attribute 
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Figure 4 | Polarisation relaxation in the SrTa0 2 N epitaxial thin films, (a) Topographic and (b) piezoresponse (A cos0) images obtained 
simultaneously from the SrTa0 2 N thin film (thickness = 288 nm) in the virgin state, (c) Phase images of the SrTa0 2 N film just after the poling process, 
(d) Averaged line profile of the dashed rectangle region in (c), which displays a clear 180° inversion of the phase, (e)-(i) Piezoresponse images obtained 
from the same area as (c) 1, 2, 4, and 8 hours after the poling process. 



the classical ferroelectric behaviour in the small domains to trans- 
type N ordering. 

Next, we validated the epitaxial stabilisation of the polar trans-type 
structure by using DFT calculations, as suggested by Page et al 6 but 
not yet confirmed. We compared the energies of the three represent- 
ative phases suggested by Page et al: paraelectric antipolar cis-type 
(space group symmetry: Pbmm), nonpolar trans-type (I4/mcm), and 
ferroelectric polar trans-type (P4mm) configurations (Table l) 6 . 
Among these structures, the cis-type Pbmm phase was the most 
energetically stable, consistent with Page et al. 6 , although we could 
not completely reproduce the absolute values of the energy difference 
among these structures (see Supplementary Note). Figure 7 shows 
how the epitaxial stress affected the phase stability, calculating the 
potential energies as functions of the in-plane lattice constant 
(Supplementary Fig. S5). Remarkably, as the in-plane lattice para- 
meters decreased, the paraelectric cis-type Pbmm phase became 
much less stable, regardless of the crystallographic orientation, while 
the stability of the trans-type structures only decreased a small 
amount. This finding implies that the compressive strain along the 
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Figure 5 | Temperature dependence of the dielectric constant of the 
SrTa0 2 N epitaxial thin film. The dielectric constants and tanc5 of the 
SrTa0 2 N film (thickness = 288 nm) measured at 20 Hz, 100 Hz, 1 kHz, 
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in-plane direction stabilised the polar trans-type nitrogen arrange- 
ment relative to the cis-type. 

Note that the cis-type Pbmm-X phase was more stable than the 
trans-type phases, even under compressive strain (Fig. 7). 
Considering the major distribution of the relaxor-like phase in the 
SrTa0 2 N film (Fig. 3b and 4b), the relaxor-like phase may be attrib- 
utable to the Pbmm-X phase or to a mixture of the cis-type phases; we 
attribute the minor classical ferroelectric phase to the trans-type 
P4mm phase. Contrasting the classical ferroelectricity originating 
from atomic displacement in the trans-type phase, the origin of the 
relaxor-like properties in the cis phases is still unclear. However, the 
lack of long-range atomic order, a common characteristic of relaxor 
ferroelectric materials 20 , may play a role in the relaxor-like behaviour. 
These results will stimulate future work to explain the relationship 
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Figure 6 | Temperature dependence of the lattice constants of the 
SrTa0 2 N epitaxial thin film, (a) The a- and c-axis lattice constants of the 
SrTa0 2 N film (thickness =100 nm) capped with an amorphous LaA10 3 
protection layer and (b) simultaneously measured lattice constant of Nb- 
doped SrTi0 3 substrate as a function of temperature. 
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Pbmm phase, E, is also displayed. Structural images are shown in 
Fig. 7 
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between the anion arrangement in perovskite oxynitrides and their 
ferroelectric properties. 

In conclusion, by using PFM we observed classical ferroelectric 
behaviour in compressively strained SrTa0 2 N epitaxial thin films 
grown on STO substrates. To our knowledge, this report provides 
the first observation of ferroelectric behaviour in perovskite oxyni- 
trides. Classical ferroelectricity coincident with visible light absorp- 
tion is difficult to achieve with conventional perovskite oxides, 
making SrTa0 2 N useful for ferroelectric photocatalysts 24 and ferro- 
electric-based photovoltaic cells 25 . We attributed the classical ferroe- 
lectricity in regions with lO^-lO 2 nm domains to trans-type anion 
ordering, likely stabilised by compressive epitaxial strain. In contrast, 
the surrounding matrix region showed relaxor ferroelectric-like beha- 
viour, characterised by the absence of spontaneous polarisation, pos- 
itive de'/dT, and a ferroelectric-paraelectric phase transition around 
600 K. DFT calculations of the phase stability under epitaxial strain 
suggested that small domains of the trans-type classical ferroelectric 
phase were embedded in a ds-type relaxor-like matrix. To further 
study the mechanism of the observed ferroelectricity, identifying the 
N-site in the oxynitrides would be a challenging study. The present 




In-plane lattice parameter (A) 

Figure 7 | Phase stability ofSrTa0 2 N. The total energy per formula unit of 
each structure, E y under epitaxial strain; Pbmm cis (diamonds, the [110] 
and [001] directions were assumed as the out-of-plane direction for 
Pbmm-X and Pbmm-X, respectively), 14/ mem trans (triangles), PAmm 
trans (circles). Out-of-plane lattice constants were optimised for each 
point and the values of E are relative to the most stable Pbmm phase 
without strain. The structures for each model without strain drawn by 
VESTA 27 are also displayed. 



study paves the way towards tailoring the functionality of perovskite 
oxynitrides by using epitaxial strain to modify anion arrangements. 

Methods 

Sample preparation. SrTa0 3 _ x N x thin films were prepared on the (100) plane of 
SrTi0 3 single crystalline substrates by NPA-PLD using a KrF excimer laser. A Sr 2 Ta 2 0 7 
ceramic pellet was used as a target, which was prepared by a conventional solid-state 
reaction of Ta 2 0 5 and SrC0 3 26 . Nitrogen gas was activated by an electron cyclotron 
resonance (ECR) plasma source (Tectra, Gen2). The partial pressure of nitrogen gas 
(P N2 ) was varied from -1.00 X 10" 7 to 1.00 X 10" 4 Torr, and the substrate 
temperature was set in a range from 500 to 800°C. The film growth rate was controlled 
by the laser repetition rate (0.3-20 Hz) and pulse energy (4-20 mj). In the PFM and 
dielectric capacitance measurements, SrTa0 2 N films deposited on conductive Nb(0.5 
wt%)-doped SrTi0 3 substrates were used after the following two-step heat treatments 
for suppressing the leakage current: (1) annealing at 500°C for 30 h under a gas mixture 
of nitrogen (P N2 = 1.0 X 10" 5 Torr, activated by ECR) and oxygen (partial pressure of 
2.0 X 10" 5 Torr) and (2) annealing at 320°C for —500 h in air. The resultant changes in 
nitrogen content and lattice parameters produced by these treatments were negligible 
within the detection limits of NRA and XRD, respectively. To avoid oxidation, we 
capped the SrTa0 2 N epitaxial thin films used for high-temperature XRD with an 
amorphous LaA10 3 layer (thickness of —20 nm). The cap layer was deposited by PLD 
at room temperature under a partial oxygen pressure of 3.1 X 10" 3 Torr. 

Chemical composition analysis. The nitrogen content in the films was evaluated by 
the NRA method using the 15 N(p,ay) 12 C resonant nuclear reaction at 898 keV. NRA 
measurements were carried out with a 1-MV tandetron accelerator at Tandem 
Accelerator Complex, University of Tsukuba. The emitted y-rays were detected with 
two 3-inch Bi 4 Ge 3 0 12 scintillation counters. The total amount of nitrogen was 
calculated by comparison to a reference sample (TiN thin film) under the reasonable 
assumption that the 15 N/ 14 N ratio was equal to the natural isotope ratio in both the 
samples and the reference. The stoichiometry of the Sr and Ta in the samples, 
fabricated under optimised conditions, was confirmed by RBS. 

Crystal structure analysis. The lattice parameters of the thin films were evaluated by 
XRD measurements with a four-axis diffractometer (Bruker AXS, d8 discover). A 
domed hot stage (Anton Paar, DHS 1100) evacuated with a rotary pump (—10 Pa) 
was used in the high temperature XRD measurement. A 200-kV TEM (JEOL Ltd., 
JEM-2010HQ and an aberration-corrected STEM (JEOL Ltd., JEM-2100F) were 
used for cross-sectional observations of thin specimens prepared by ion milling. 

Optical absorption spectrum. The absorption spectrum of the SrTa0 2 N films was 
calculated from extinction coefficients evaluated by a spectroscopic ellipsometer (J.A. 
Woollam, M-2000U). 

Atomic force microscopy (AFM) and PFM. AFM and PFM measurements were 
conducted using a scanning probe microscope (SH-nanotechnology, E-Sweep with 
NanoNaviReal) with a homemade PFM system composed of a function generator 
(NF corporation, 1930A) and a lock- in amplifier (EG&G, 7260). An Rh-coated silicon 
cantilever (SH-nano technology, SI-DF3-R(100) with stiffness of 1.1 N/m and a 
resonance frequency of 25 kHz in air) was employed for the measurements. In the 
PFM measurements, the frequency of the driving AC electric field was set to a value 
slightly (less than 1 KHz) deviated from the contact resonance (—110 kHz). The 
amplitude of the driving AC field was 1 V. A domain poling process was performed 
using the following three steps: (1) applying a sample bias of +7 V on the central 
5X5 urn 2 area, (2) applying —7 V on the central 3X3 um 2 area, and (3) applying 
+ 7 V on the central 1X1 urn 2 area. 

Dielectric measurements. Dielectric constants were determined by a conventional 
capacitance method. Pt top electrodes (200 \im(p) were sputtered onto the SrTa0 2 N 
thin films. The capacitance and tanc) of the SrTa0 2 N films were examined with a 
precision LCR meter (Agilent, E4980). A cryogen-free micro -manipulated probe 
station (Lake Shore, CRX-4K) was used for temperature control. 

DFT calculations. The potential energy surface (PES) was calculated using Vienna ab 
initio software package. The Perdew-Burke-Ernzerhof (PBE) functional and 
projector-augmented wavefunction pseudopotential were adapted. The cut-off 
energy was 560 eV (41.2 Ry). The k-points were sampled according to the 
Monkhorst-Pack (MP) 8 X 8 X 8 for Pbmm and PAmm cells. The MP 8 X 8 X 4 was 
used for 141 mem. The atomic positions were optimised until the residual force acting 
on all atomic coordinates was less than 0.01 eV/A. The bottom of the PES was 
calculated by optimising the size and shape of the unit cell and the atomic positions. 
The residual pressure was smaller than 0.035 GPa in all directions for all considered 
anion orderings. A validity check was performed using the Quantum Espresso (QE) 
software package. In the case of QE, the PBE functional and ultrasoft pseudopotential 
were employed, and the cut-off energy was 48.0 Ry for the wavefunctions. The k- 
points were sampled according to MP 4X4X6, MP 6X6X6, and MP 6 X 6 X 4 for 
Pbmm, P4mm, and 14/mcm, respectively. 
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